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The resistivity as function of temperature of high temperature superconductors is very unusual and despite its 
importance lacks an unified theoretical explanation. It is linear with the temperature for overdoped compounds 
but it falls more quickly as the doping level decreases. The resistivity of underdoped cuprates increases as an 
insulator below a characteristic temperature where it shows a minimum. We show that this overall behavior 
can be explained by calculations using an electronic phase segregation into two main component phases with 
low and high electronic densities. The total resistence is calculated by the various contributions through several 
random picking processes of the local resistivities and using a common statistical Random Resistor Network 
approach. 

PACS numbers: 74.20.-z, 74.25.Dw, 74.72.Hs, 64.60.Cn 



I. INTRODUCTION 

The search for the underlying mechanism of the High crit- 
ical temperature superconductors (HTSC) has become one of 
the most studied problems in condensed matter physics. De- 
spite the intense experimental and theoretical effort there is 
neither a consensus on the basic mechanism nor on the the 
generic phase diagram. In opposition to the low temperature 
superconductors that the normal phase is a well known Fermi 
liquid, the normal phase of HTSC has many non understood 
and non Fermi liquid like properties. For instance, they have 
a normal state gap or pseudogap that in some cases remain at 
temperatures much higher than the resistivity transition tem- 
perature Tc{l] and whose origin is still a matter of great de- 
bate. 

The pseudogap manifests itself by a suppression of spectral 
weight of the normal- state electronic density of state as es- 
tablished by many different experiments on bulk and surface, 
momentum and real space, as has been extensively discussed 
in many review articles[l- 3]. Consequently, it has become ac- 
cepted that the HTSC phase diagram is dominated by two en- 
ergy scale and two distinct quasiparticle dynamics yl 3]. One 
associated with the superconducting phase the other associate 
with the pseudogap temperature T*. 

Another important feature of the cuprate superconductors 
is the presence of inhomogeneities. A phase separation tran- 
sition was observed experimentally a long time ago in the 
La2Cu0^j^y^. In this experiment the phase separation of 
the interstitials oxygen are accompanied by an electronic 
phase separation of holes in CUO2 planes. Many different 
experiments have added evidences that the charge distribu- 
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tion in the CUO2 planes of the high temperature superconduc- 
tors (HTSC) is microscopically inhomogeneous like neutron 
diffractionfl |7|], muon spin relaxation (/i5'7?)[8], NQR and 
NMrI^ [i3] . These experiments indicated that the inhomo- 
geneities are more pronounced on the underdoped side of the 
phase diagram and related it with the non Fermi liquid behav- 
ior of the normal phase. The subject of phase separ Recent 
scanning tunneling microscopy (STM) studies on Bi2212 re- 
veal spatial variations of the electronic gap amplitude on a 
nanometer length scale even on overdoped compounds flTI- 
Il3l1 . The local densit y of states (LDOS) measured in these 
STM experiments |[ll[4l3|] has two different forms that are 
possibly associated with a segregation into metallic and in- 
sulator nanoscopic regions |[l4ll . Consequently the electronic 
phase separation in cuprates has been the subject of many 
articlesilsl-lIS] and many books iBB. 

The electronic disorder of the normal phase are manifested 
in many transport properties of a HTSC series, like the re- 
sistivity as function of the temperature R{T), that despite its 
importance, lacks a widely accepted explanation. Here we 
want to show that the particular behavior of R{T) for dif- 
ferent compounds of a single series is a manifestation of 
the intrinsic inhomogeneity of these materials. We can di- 
vide the resistivity behavior for a system in three general 
types: It falls linearly with the temperature for samples in 
the overdoped regions i26l] . It falls below the linear trend for 
compounds near the optimum values as measured by several 
groups i!?! |28t]. For very low doping it falls down with the 
temperature, reaches a minimum at an intermediate temper- 
ature, starts to increase at lower temperature until it reaches 
a local maximum, an then it drops to zero below Tc. This 
such a reentrant behavior is observed in almost all underdoped 
cuprates i29l - [3lll but its origin has not been very much dis- 
cussed in the literature. The significance and origin of such 
rich phenomena will be the focus of our work. 

Based on the above facts we perform calculations taking 



into account the EPS transition by the use of the general phase 
separation theory of Cahn-Hilliard (CH) |[32ll . It describes how 
a system evolves from small fluctuations around an average 
density /? to a complete separation into low and high density 
regions. At each of these small regions we will associated 
a local resistivity derived from the work of Takagi et al ll26ll . 
who performed systematic measurements on the LSCO series 
(La2_pSrpCu04) to a large range of values of T and p. Then 
we use the Random Resistor Network (RRN) methodi33] to 
derive the system resistivity as a function of the temperature 
R{T) in a similar way as was applied to the manganites[34]. 
The local density and resistivity is randomly picked and the fi- 
nal 7? (r) result is an average of many different configurations. 
We show that the non conventional features discussed in the 
former paragraph are well described by this approach. 



II. ELECTRONIC PHASE SEPARATION 

The electronic phase separation was found by several 
HTSC experiments and detected to be in the form of either 
stripes 10, [35I], patchwork ll36|] or checkerboard fs?!]. Earlier it 
was believed that the underline mechanism of electronic phase 
separation was associated with defects, oxygen or cation dis- 
order. Indeed, oxygen phase segregation was observed on 
the La2Cu04^^ by x-ray and transport measurements lisi flsll . 
However depending on the type of experiment, there are 
some HTSC materials, specially YBCO, that appear to be 
more homogeneous or, at least do not display any gross 

inhomogeneity (M El] • 

On the other hand, electronic phase separation appear to be 
an universal feature of HTSC. It is seen in the entire series of 
LSCO by neutron diffraction JtI] and in many BSCO samples 
by STMlfrTI-fBll. ARPES have also detected distinct quasipar- 
ticle behavior in LSCO and BSCO. Recently inhomogeneous 
magnetic-field response measurements on some LSCO com- 
pounds and on the YBa2Cu30y (YBCO) family ||4Q|] were also 
interpreted as compatible with electronic phase separation in 
form of domains. Specially in the case of YBCO where dop- 
ing occurs via a change in the oxygen concentration of the 
CuO-chain, the inhomogeneities cannot be attributed to dis- 
order in the cations substitution, providing strong support to 
an intrinsic EPS transition. The effect of the annealing time 
connected to the structural organization of oxygen intersti- 
tials in La^CuO/^-y was studied recently X-rays diffraction 
techniques 1I22I] 

The possible origin of this EPS transition is the proxim- 
ity to the insulator AF phase, common to all cuprates, as we 
derived from the principle of the competing minimum free 
energy IH]. These calculations yield a line Tps{p) close to 
the experimental results reported by Timusk and StattiH] (T^ 
in their notation), and it has been taken as a linear function 
just for simplification as it is schematically shown in Fig.© 
and dU. When the temperature decreases below Tps{p), the 
free energy of the homogeneous system with average density 
p becomes lower than the anisotropic one made of a bimodal 
distribution |[T6|] of AF domains with p{i) ^0 and high hole 
density domains with p{i) ^ 2p. This approach is justified 



by the NQR experimental dataf^] which was interpreted as 
coming from two distinct regions with symmetric local dop- 
ing around the average value p. The two symmetrical local 
dopings increase their differences as the temperature is low- 
ered, that is characteristic of second order transitions. 

To describe the formation of the local hole density domains 
we define the difference between the local and the average 
charge density u{i^T) = {p{i^T) — p) as the EPS transition 
order parameter. Clearly u{i^T) = corresponds to the ho- 
mogeneous system above Tps{p). Then the typical Ginzburg- 
Landau free energy functional in terms of such order parame- 
ter near the transition is given by 

= y^\Vu{iJ)\^^V{u{i,T)). (1) 

Where the potential V{u,T) = A^{T)u^ /2 ^ B^u^ /A ^ 
A^{T) = a{Tps{p) — r), a and B are constants that lead to 
lines of constant values of A{T)/B, parallel to Tps{p), as 
shown Fig([T]). £ gives the size of the grain boundaries among 
two low and high density phases ;?±(/)|[i7l Ull]. The en- 
ergy barrier between two grains of distinct phases is Eg (T) = 
A^{T) /B that is proportional to {Tps — T)^ near the transition, 
and becomes nearly constant at low temperatures. The two 
minima of the y(w, T) is shown in Fig.([Tl 
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FIG. 1 . The potential commonly used in the free energy of the CH 
equation which gives rise to phase separation as function of the order 
parameter u or the density p. Notice that the two minima yield the 
two equilibrium low and high density phases and the energy barrier 
or free energy wall is Eg and depends on how low the temperature is 
below Tps, that is, the temperature difference Tpg — T. 

Thus, hereafter we will use Eg{p^T) =V{p^T) sls the grain 
boundary potential or the binding potential for holes in a given 
grain. V{p^T) =V{p) xV{T) and we assume, for simplicity, 
that V{p) has a linear behavior on p, whose equipotentials are 
parallel to Tps{p). In Fig.© we plot Tps{p) where the phase 
separation starts, To{p) where the EPS is detected and Ti{P) 
where insulators grains appears, all assumed to be linear on 
p just for simplification but in agreement with some signal 
detected by many experiments DM and DI are phases 



with low disorder in the sense that the two phases are metalHc 
(DM) or insulator (DI). At lower temperatures, the disorder 
increases and the system is composed of very low density or 
insulator regions and very high or metallic regions (DM+DI). 
See Fig.(|2l 

The CH equation can be written 1I42I1 in the form of a con- 
tinuity equation of the local density of free energy /, dfU = 
— V • J, with the current J = MV(5// 5u), where M is the mo- 
bility or the charge transport coefficient. Therefore, 



^ = -MV\e^V^u^A\T)u-B^u'). 
ot 



(2) 



In previous pap ers we have verified the applicability of the 
CH approach 14 ill . We have made a detailed study of the den- 
sity profile evolution in a 105 x 105 array as function of the 
time steps, up to the stabilization of the local densities, for pa- 
rameters that yield stripe ifTsIl and patch work|[l7[ E] patterns. 

Here we want to concentrate on the energy structure of the 
grains which leads to additional scattering and changes the 
whole charge dynamics. We show that below the Tps, the 
grain boundary energies Eg (see the inset of Fig.Q) increase 
as the temperature goes down, segregating and confining holes 
at certain regions or grains with an increasing difference in the 
local doping. The confining potential favors pair formation at 
low temperatures, giving rise to the local superconductivity in 
isolated regions. 
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FIG. 2. Important crossover lines derived from many experiments (il] 
that are associated with the degree of phase separation. The elec- 
tronic segregation starts at Tps{p) resulting into two types of phases, 
high and low densities as a consequence of the potential shown in 
Fig.l. DI and DM stands for disordered insulator and disordered 
metal phases. DI+DM stands for a larger disorder with insulator and 
metallic phases cohexisting. DMi -\-DM2 is the situation which both 
phases are metallic, but with different densities. 

As the temperature goes down below the Tps{p), the sys- 
tem increases the phase separation and in favorable conditions 
tends to form a bimodal distribution made of low and high 
density domains. Fig.© shows the various levels of phase 
separation as function of the temperature difference | Tps — T\. 
At low density (/? < 0. 16) and low temperature the system can 
achieve the maximum charge segregation and become com- 




posed essentially of two densities, one of p{i) ^ and other 
of p{i) ^ 2p. There are disordered insulator (DI) regions and 
disordered metallic (DM) ones. In Fig.© we show this situ- 
ation for the case with /? = 0.16 and T ^ OK, with half of the 
sites with p{i) ^0 and the other half with p{i) 0.32, as it is 
seen in the inset histogram made over all the sites " / " of the 
system. 

The density map shown in Fig.© is made on a 100 x 100 
unit cells square lattice, the low and high density grains have 
different size and are composed by 30-100 cells. The inset is 
the histogram of the local doping at each crystalline site that 
shows how the density segregates into two main bands whose 
width values depend on temperature. 

At low temperatures below To{p) (see Fig.©), we can 
consider the grains as isolated entities with very weak inter- 
grain coupling and a typical Fermi energy of a grain is much 
smaller than the bulk Debye energy, the so called anti adia- 
batic regime [44]. 

Thus below the Tq line in Fig.©, a cuprate is composed of a 
mixture of insulator (p{i) ^ 0) and metallic (p{i) ^ 2p) grains. 
A clear experimental realization of such system is the granular 
Pb film on a glass substrate which can change from insulator 
to superconducting behavior with increasing Pb grains ll45ll . 
For very low average doping the system is mainly insulator, 
but for > 0.05 the large density grains become metallic like, 
and in general both types of regions are present in one single 
sample. 




FIG. 3. (color online) Density map of system at low temperatures. 
The low and high density domains are clearly formed. The inset is 
the local density p{i) histogram with average hole density p = 0.16, 
showing that the high density grains are metallic like and the low 
density (insulator) ones are made of p{i) ^ 0. 

As the temperature decreases below Ton{p) (see Fig|4]) it 
is possible to show that the d-wave superconducting ampli- 
tude Aj(/, r) appears in some isolated metallic grains lfl9ll . 
what is known as intragrain superconductivity by the the- 
ory of granular superconductors ||45l. And the zero resistivity 
transition takes place only when the Josephson coupling Ej 
among these grains is sufficiently large to overcome thermal 
fluctuations 1 45], that is, Ej{p^T = Tc) ~ kBTc{p) what leads 
to phase locking and long range phase coherence. 



In the next section we will use the above results to calculate 
the average resistivity. The main point is that any compound 
of a HTSC system goes through a EPS transitions and the lo- 
cal densities follow a bimodal distribution with a dispersion 
similar as that shown in the insert histogram of Fig.®. This 
is interpreted as a distribution of local resistivities and will be 
used in connection with the RRNIiSSi]. 

III. THE RESISTIVITY CALCULATIONS 

The RRN has been quite useful in the study of systems with 
this type of disorder like on the percolation phenomenon as 
well as conductance models themselves Ji^. Here we adopt 
a square lattice, where nodes are connected by resistive links, 
which can have different resistances. Each of these would 
be a nanoscopic domain inside our sample, which could be 
metallic, insulating (AF phase) or superconducting. The dis- 
tribution of their fractions depend on temperature and nominal 
doping composition. Two bus bars are imagined to be placed 
at top and bottom of the array so that a external voltage source 
can set a voltage drop V across the network. The total resis- 
tivity will be simply obtained by the ratio V/I where / is the 
total current through circuit. 

As we discussed in the previous section, above the phase 
separation line (PS), the system has only metallic links, and 
resistance is considered to be a linear function of temperature. 
Below PS, phase separation begins and two bands of compo- 
sition p are allowed. These bands have their central values Ap 
apart and are 5p{i) wide. Both Ap{i) and 5p{i) depend on 
temperature and are all symmetric to the nominal or average 
composition p. In Fig.Q the bands are drawn in a schematic 
phase diagram for a nominal doping fraction p = 0.01. The 
temperature evolution of these bands follows the study of the 
histograms at different temperatures as the one shown in the 
inset of Fig.©. 
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FIG. 4. (color online) Schematic phase diagram showing the EPS 
with bands of doping fractions Sp(i) as a function of temperature. 
At low temperature below Ton tiny domains of SC appear on the right 
band, and for low doping compounds small domains of AF insulators 
appear on the left band inside the sample. 

Resistivity is a function of both composition and tempera- 



ture. From the results of Takagi et al |[26h . for various com- 
positions p and a given temperature, is possible to devise a 
nearly exponential dependence on composition. The metallic 
links in the RRN will then follow a random distribution in the 
allowed values of p{i) (the bands) and their resistances will be 
given by a function derived from the values of Takagi et al ll26l] 

R{p{i)J) =A{T)cxp{-B{p{i)-p)} (3) 

Where B(T)=0.05 and A(T) are derived directly from the 
LSCO series measurements ll26ll . 

On the other hand, following experimental results Hsl |28l 
45] and calculations [3 similar to the intragrain super- 
conductivity discussed in the previous sections, some super- 
conducting regions appear at low temperature. Thus, as the 
temperature decreases below Ton{p), a fraction of the most 
conducting metallic links (largest ;?(/)'s) are replaced by SC 
links and this fraction reaches 100% of the metallic band at 
Tc{p). Similarly, insulators appear in the low density band as 
p^O, which happens at lower temperatures, as can be seen in 
Fig.©. Likewise the CH phase separation calculations shown 
in the previous section, the resistivity calculations were made 
on 100 X 100 lattices. Periodic boundary conditions were used 
only on the sides without the V source terminals. The result- 
ing linear systems that come from Kirchhoff laws were solved 
with the aid of the subroutine ina57, provided by HSL, that 
makes use of a multifrontal algorithm JItIi, what speeded up 
the calculations. To each value of p and T we use 300 samples 
in order to obtain results independent of the random process 
of sharing the p{i) distribution. 



— Naqib, 2005 
o simulation RRN 

1,5- 




100 200 300 



T(K) 

FIG. 5. (color online) Simulation results compared with experimen- 
tal data for the case Yq ggCaQ 2oB^2^%^7-5 ^^^^ average doping 
p = 0.136. As the temperature decreases below Ton = \60K some of 
the resistance become superconducting and total resistivity decreases 
below the linear regime. An average with 300 samples were used for 
each point in the graph. 

Results of simulations for typical samples near op- 
timum doping are shown in Fig|5] for the compound 
Yq ggCaQ 2oBa2Cu30^_^ with nominal doping p = 0.136. 
300 samples were used for each point in the graph. The 
phase separation begins at 400 K; SC links of resistivity 



1 X 10~^m^l-cm are introduced below 180K, according to 
the Ton{p) line. 

The small error bars can tell how reliable the results are, 
that is, a very small dispersion is attained with averages made 
on 300 configurations. It is also remarkable how the p x T 
linear dependence is attained in the calculations with the two 
bands scenario. That can be explained by the proper choice 
of the function R{p^ T) defined in Eq. (O, which can be bet- 
ter understood under the light of results obtained by Costa et 
ah who used renormalization group theory to study square 
When a square RRN has resistors of only two con- 
ductivity values g\ and g2, distributed with probabilities q—\ 
and q, respectively, Costa et ah showed that the network over- 
all conductance, expressed as a function of probability q, is 



a{q) a{l-q) ^ gi 
CT(1) g2 



(4) 



With equal probabilities, i.e. ^ = 0.5, its overall resistance 
will simply be 

r y/nri. (5) 

In our case, if the widths bp were chosen to be zero, a 
network of links with resistances R{jp — lS.p{i)^T) and R{p + 
Ap{i)^T) would result in an overall resistance R{T) = A{T), 
easily derived from Eqs. ^ and (O. The same result holds 
for finite width 5 p{i), if symmetric distributions around p are 
used. Making A (r) linear is then the right choice. Departures 
from this linear behavior are explained by the growth of SC 
and AF domains. 

The exact linear dependence of A(r) for a given doping 
fraction p can be extracted from the linear part of experimen- 
tal curves as in Fig s. ([5] and O 

The fraction of SC links to be inserted below Ton was also 
calculated from Eq. IH as follows. The difference between 
function A (r) and experimental data was normalized by A(r) 
itself, and Eq. |4] gave the probability q associated with a 
network of conductances gi = I and g2 = I x 10+^m^l-cm. 
When distributing p{i) in the allowed bands by means of ran- 
dom numbers x, uniformly distributed in the interval [0,1], 
those links for which x> {I — q) (the best conductors) were 
made superconductors. 

The resistivity for low doping materials requires special 
care because it may contain insulator, metallic and supercon- 
ducting regions represented by links in the RRN method. Thus 
simulations for a case which reentrant behavior is present 
were made for the p = 0.07 compound and the results are 
shown in Fig. [6l We find that such complex behavior is due 
to the mixture of superconducting and insulator regions in the 
two different bands. The appearance of the insulator regions 
in our RRN matrix occurs when 5p{i) crosses the value of 
p{i) < 0.05 JsUl at low temperatures as it is demonstrated in 
Fig.[6l This is the reason why, at low temperatures, the resis- 
tivity increases, but eventually, as the temperature decreases 
further and crosses Ton, the metallic resistors become super- 
conductors and the resistivity vanishes due to the percolation 
among the superconducting regions, as we explain below. 

The appearance of the superconducting links (resistivity 
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FIG. 6. Simulation results (circles) for in-plane resistivity of 
La2_pSrpCu04 with p = 0.07 and the reentrant behavior. Each cal- 
culated point is the average of 600 configurations. Dots represent 
experimental data obtained by Takagi et al. 123] . 



1 X 10~^ mO. • cm) begins when the temperature decreases be- 
low 280 K. Their fraction is derived in a similar fashion, 
based on Eq. Q, as described above. Insulating links (p = 
200m^l-cm) were first added at 140 K as p{i) reaches the 
value 0.05 in the low density band. In this case, due to the 
low values of the densities, the number of samples for each 
point was 600 for temperatures below 60 K and 300 for the 
remaining temperatures. As T ^ both insulating and SC 
links fractions come close to the link percolation threshold 
0.50. That is the reason for the strong variation expressed in 
the error bars as T ^ in Fig. [6l 



IV. CONCLUSION 

In this paper we derived the resistivity of a complete HTSC 
series assuming an EPS transition. Following many experi- 
mental results we took into consideration that a HTSC com- 
pound undergo an electronic segregation and is composed 
mainly of two types of regions or grains. These are the high 
and low local density which gives a metallic or insulator be- 
havior to the grains as observed by the STM measurements of 
the LDOS. We used simulations with the CH theory to follow 
the effects of the EPS as function of the temperature in order 
to describe the disordered composition of a HTSC sample. 

With this approach the overall resistivity can be calculated 
through the RRN method [33] and used to reproduce the mea- 
sured values either in the overdoped or underdoped regions. 
Our method based on the effects of disorder is in very good 
agreement with the measured departure of the linear behav- 
ior of near optimum compounds and the completely differ- 
ent down and up reentrant behavior for weakly doped samples 
that have been studied in detail! 31], but so far without a sim- 
ple and unified interpretation. 

The excellent fitting on such a rich variety of features with 
one single approach without adjustable parameters and with 
values of the local resistivities taken from the measurements 
of Takagi et al ll26|] for an entire series, let us conclude the 



following: the EPS described here and the onset of local orx 
intragrain superconductivity given by the curve Ton are univer- 
sal properties and must be considered to interpret the transport 
properties of any HTSC system. 



V. ACKNOWLEDGMENT 



We gratefully acknowledge partial financial aid from 
Brazilian agencies CNPq and Capes. 



[1] T. Timusk and B. Statt, Rep. Prog. Phys., 62, 61 (1999). 

[2] J.L. Tallon and J.W. Loram, Physica C 349, 53 (2001). 

[3] S. Huffner, M.A. Hossain, A. Damascelli, G.A. Sawatzky, Rep. 

Prog. Phys. 71, 062501 (2008). 
[4] M. Le Tacon, A. Sacuto, A. Georges, G. KotHar, Y. Gallais, D. 

Colson and A. Forget, Nature Physics 2, 537 (2006). 
[5] J. D. Jorgensen, B. Dabrowski, Shiyou Pei, D. G. Hinks, L. 

Soderholm, B. Morosin, J. E. Schirber, E. L. Venturini, and D. 

S. Ginley Phys. Rev. B38, 11337 (1988). 
[6] J.M.Tranquada, B.J. Sternlieb, J.D. Axe, Y. Nakamura, and S. 

Uchida, Nature (London),375, 561 (1995). 
[7] E.S. Bozin, G.H. Kwei, H. Takagi, and S.J.L. Billinge, Phys. 

Rev. Lett. 84, 5856, (2000). 
[8] Y. J. Uemura, Sol. St. Phys., 126, 23 (2003). 
[9] P M. Singer, A. W. Hunt, and T. Imai, Phys. Rev Lett. 88, 

47602 (2002). 

[10] H.-J. Grafe, N. J. Curro, M. Hucker, and B. Biichner, Phys. Rev 

Lett., 96, 017002 (2006). 
[11] K. McElroy, D.-H. Lee, J. E. Hoffman, K. M Lang, E. W. Hud- 

son, H. Eisaki, S. Uchida, J. Lee, J.C. Davis, c ond-mat/04040051 

and Phys. Rev Lett. 94, 197005 (2005). 
[12] Kenjiro K. Gomes, Abhay N. Pasupathy, Aakash Pushp, Shim- 

pei Ono, Yoichi Ando, and Ali Yazdani, Nature 447, 569 

(2007). 

[13] Abhay N. Pasupathy, Kenjiro K. Gomes, Colin V.Parker, Jin- 
sheng Wen, Zhijun Xu, Genda Gu, Shimpei Ono, Yoichi Ando, 
Ah Yazdani, Science 320, 196 (2008). 

[14] E. V. L. de Mello, R.B. Kasal, J. Supercond. Nov Magn., 24, 
1123 (2011). 

[15] J.C. Grenier, N. Lagueyte, A. Wattiaux, J.P. Doumerc, P. Dor- 
dor, J. Etourneau, M. Puchard, J. B. Goodenough, and J.S. 
Zhou, Phys. C202, 209 (1992). 

[16] E.V.L. de Mello, E.S. Caixeiro, and J.L. Gonzalez, Phys. Rev. 
B67, 024502 (2003). 

[17] E.V.L. de Mello, and E.S. Caixeiro, Phys. Rev B70, 224517 
(2004). 

[18] E. V. L. de Mello, and D. N. Dias, J. Phys. CM. 19, 086218 
(2007). 

[19] E. V. L. de Mello, R. B. Kasal, C. A. C. Passos, J. Phys,: Con- 
dens. Matter 21, 235701 (2009). 

[20] D. Innocent! et al, J. Supercond. Nov Magn., 22, 529 (2009). 

[21] K.I. Kugel et al, Phys. Rev B78, 165124 (2008). 

[22] M. Fratini, N. Poccia, A Ricci, G. Campi, M Burghammer, G. 
Aepph, and A. Bianconi, Nature 466, 841, (2010). 

[23] "Phase Separation in Cuprate Superconductors", edited by E. 
Sigmund and K.A. Muller (Springer- Verlag, Berlin, 1994). 



[24] "Stripes and Related Phenomena", edited by A. Bianconi and 
N. L. Saini, (Academic/Plenum publisher. New York, 2000). 

[25] "Symmetry and Heterogenity in high temperature Supercon- 
ductors", edited by A. Bianconi. Nato Science Series II 
(Springer Dordrecht, Netherlands, 2006) 

[26] H. Takagi et al , Phys. Rev Lett. 69, 2975 (1992). 

[27] S. H. Naqib et al, Phys. Rev B71, 054502 (2005). 

[28] C. A. Passos, M. T. Orlando, J. L. Passamai, E. V. de Mello, H. 
P Correa, and L. G. Martinez, Phys. Rev B 74, 094514 (2006). 

[29] Y. Ando et al, Phys. Rev Lett. 75, 4662 (1995). 

[30] S. Ono et al , Phys. Rev Lett. 85, 638 (2000). 

[31] Seongshik Oh et al, Phys. Rev. Lett. 96, 107003 (2006). 

[32] J.W. Cahn and J.E. Hilhard, J. Chem. Phys, 28, 258 (1958). 

[33] S. Kirkpatrick, Rev Mod. Phys. 45, 574 (1973). 

[34] Matthias Mayr,Adriana Moreo, Jose A. Vergilio, Jeanette 
Arispe, Adrian Feiguin, and Elbio Dagotto, Phys. Rev. Lett. 86, 
135 (2001) 

[35] A. Bianconi, N.L. Saini, A. Lanzara, M. Missori, T. Rossetti, 
H. Oyanagi, H. Yamaguchi, K. Oda and T. Ito, Phys. Rev. Lett. 
76, 3412 (1996). 

[36] S. H. Pan et al.. Nature, 413, 282-285 (2001) and 

cond-mat/0107347 
[37] T. Hanaguri, C. Lupien, Y. Kohsaka, D. -H. Lee, M. Azuma, M. 

Takano, H. Takagi, J. C. Davis, Nature 430, 1001 (2004). 
[38] J. Bobroff, H. Alloul, S. Ouazi, P Mendels, A. Mahajan, N. 

Blanchard, G. Collin, V. Guillen, and J.-F. Marucco, Phys. Rev. 

Lett.89 157002 (2002). 
[39] J.W. Loram, J.L. Tallon, and W. Y. Liang, Phys. Rev B69, 

060502 (2004). 
[40] J. E. Sonier, et al, Phs. Rev Lett. 101, 117001 (2008). 
[41] E.V.L de Mello, and Otton T. Silveira Filho Phyica A 347, 429 

(2005). 

[42] A.J. Bray, Adv Phys. 43, 347 (1994). 

[43] E.V.L. de Mello, R. B. Kasal, C. A. C. Passos, Otton T. S. Filho, 

Physica B 404 3 119 (2009). 
[44] E. V. L. de Mello and J. Ranninger, Phys. Rev B55, 14872 

(1997). 

[45] L. Merchant, J. Ostrick, R. P. Barber Jr, and R. C. Dynes, Phys. 

Rev B63, 134508 (2001). 
[46] J. S. Andrade, Jr., N. Ito, and Y. Shibusa, Phys. Rev B 54, 3910 

(1996). 

[47] I. S. Duff, ACM Trans. Math. Softw. 30, 118-144 (2004). 
[48] U. M. S. Costa, C. Tsalhs, and G. Schwachheim, Phys. Rev B 
33,510(1986). 



